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ABSTRACT 

We present a radial velocity study of the RY Tau jet-disk system, designed to determine whether 
a transfer of angular momentum from disk to jet can be observed. Many recent studies report on 
the rotation of T Tauri disks, and on what may be a signature of T Tauri jet rotation. However, 
due to observational difficulties, few studies report on both disk and jet within the same system 
to establish if the senses of rotation match and hence can be interpreted as a transfer of angular 
momentum. We report a clear signature of Keplerian rotation in the RY Tau disk, based on Plateau 
de Bure observations. We also report on the transverse radial velocity profile of the RY Tau jet 
close to the star. We identify two distinct profile shapes: a v-shape which appears near jet shock 
positions, and a flat profile which appears between shocks. We do not detect a rotation signature 
above 3 sigma uncertainty in any of our transverse cuts of the jet. Nevertheless, if the jet is currently 
in steady-state, the errors themselves provide a valuable upper limit on the jet toroidal velocity of 10 
km s“^, implying a launch radius of < 0.45 AU. However, possible contamination of jet kinematics, 
via shocks or precession, prevents any firm constraint on the jet launch point, since most of its angular 
momentum could be stored in magnetic form rather than in rotation of matter. 


Keywords: ISM: jets and outflows - stars: individual (RY 
Tau) - submillimeter: ISM - infrared: ISM - ISM: kine¬ 
matics and dynamics 

1. INTRODUCTION 

The origin of protostellar jets, their role in angular 
momentum extraction and their impact on protoplan¬ 
etary disks are still critical open issues in star forma- 
tio n. The MHD ejection m echanism originally proposed 
by (jBlandford et al.lll982[) is currently the favoured pro¬ 
cess for jet generation as it can naturally account for 
both high jet speeds and high degrees of collimation 
close to the source. MHD jets from young accreting 
stars may hav e three possible components (see review by 
iFerreira et al.l (I2006D and references therein): a strong 
pressure-driven stellar wind (SW); sporadic wide-angle 
magnetospheric ejections (MEs) caused by the interac¬ 
tion of the stellar magnetosphere with the inner disk 
edge and a self-collimated centrifugally-driven disk wind 
either originating from the co-rotation radius (X-wind so¬ 
lution) or from a wider range of disk radii (extended disk 
winds). MEs and SWs currently appear as the best can¬ 


didates to spin-down young accreting stars below break¬ 
up. If extended disk-winds are present in young stars, 
this would imply that the inner astronomical units of 
proto-planetary disks are more strongly magnetised than 
assumed so far, which would have strong implications for 
angular momentum transfer and p lanetary formation and 
migration models (se e reviews by iTurner et al.ll201^ and 
IBaruteau et al.l[2014D . 

Over the last decade, both semi-analytical mod¬ 
elling and numerical simulations have been extensively 
used to investigate how the mass-loading, kinematics 
and collimation properties of disk winds are contribed 
but the detailed physics of the underlying accretion 
disk (magnetic field int ensity and topolog y , thermody¬ 
nami cs, resist iyity etc (ICasse et al.l (2000t iFendt etakl 
20021 2006. I^iirii et al 2007 : Murphy et al.l l20ICt 


zuvA ZUUD. iZiamii et ai. \zoo(t liviurDny et al.l I2uiut 
Sheikhnezami et al. I l2012t Isteoanoys et al.ll2014D . A di¬ 


rect comparison of jet observations with model outputs 
remains difficult because current numerical works remain 
biased towards small magnetic lev er arm solution s . proR 
ably due to numerical diffusivity iMurphv et al.l (120101) . 
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In addition, a proper treatment of non-equilibrium tem¬ 
perature and ionisation is required, which can com¬ 
pletely change the emission line profiles for the under- 
lying dynamics . So far, this has only been done by 
iTesileami et alJ (j2014fl . and only for simulations treat¬ 
ing the disk as a platform, which are not dynamically 
self-consistent. Conyersely, semi-analytical solutions re¬ 
main limited by the self-similar assumption and the lack 
of proper inner and outer boundary conditions. 

An alternatiye approach, which we adopt here, is 
that of relying on general conseryation principles of 
any steady, axisymmetric, cold magneto-centrifugal disk 
wind solution after it is launched and in ideal MHD, 
namely the conseryation of total (matter plus magnetic) 
specific energ y and angular moment um along each mag¬ 
netic surface. lAnderson et al.l (|2003r ) showed that a com¬ 
bination of these two inyariants allows us to get rid of the 
non-obseryable magnetic terms and constrain the launch¬ 
ing radius of the disk wind streamline independently of 
any specific MHD model, simply from the obser y ed jet 
rotational and poloidal yelocity. iFerreira et al.l (120061 1 
showed that a lower limit on the disk wind magnetic 
leyer arm can also be obtained in this way. 

In recent years, pushing the limits of obseryational res¬ 
olution has reyealed a number of exciting detections of 
gradients in th e radial velocity prof i le across jets from 
T Ta uri stars (iBacciotti et al.l l2002l iCoffev et al.l 12004 
2007. lWoitas et al.ll2005ll . These may be interpreted as a 
signature of jet rotation about its symmetry axis, thereby 
representing the long-awaited observational support of 
the theory that MHD jets extract angular momentum 
from star-disk systems. Under the assumption of steady 
mass-loss, deriving jet angular moment um has the po¬ 
tential to constrain the laun ching radius (IBacciotti et al.l 
120021: lAnderson et al.l l2003f l. If the observed transverse 
velocity gradients are due t o rotation, then they imply 
launching radii of 0.1-3 AU (|Ferreira et al.l[^006fi . How¬ 
ever, the possibility that we are indeed observing jet 
rotation in pre-mainsequence systems is undergoing ac¬ 
tive debate. Alternative interpretations in clude asym¬ 
metric shocking and/or jet precession (e.g . lSokeil[2?i?)M 
iCeroueira et al.l [200(11 iCorreia et al.l l200'^ . It is there¬ 
fore important to perform new checks on the rotation 
interpretation. 

First, to ensure that we probe the kinematics of the 
launching process, the jet must be observed as close as 
possible to the star where any evidence of angular mo¬ 
mentum transfer is still preserved, rather than far from 
the star by which time the internal shocks, or inter¬ 
action with the ambient medium, may have disrupted 
any rotation signature intrinsic to the launch mecha¬ 
nism. Second, we must check that rotation signatures 
are consistent at different positions along the jet (in 
various tracers and epochs). This has been checked in 
only three studies so far: an optica l study of DG Tau 
(|Bacciotti et al.l[2002ll and RW Aur (iWoitas et al.ll2005ll 
where consistency was found; a UV study of DG Tau 
and RW Aur where co nsistency was found in the for¬ 
mer but not the latter (jCoffev et al.l 120071 Coffe^^^jl 
2012 ); and a near infrared study of DG Tau (|White et al. 
20141 where systematic transverse velocity gradients 
could not be identified). Third, we must check that the 
jet Doppler gradient matches the direction of the disk 
Doppler gradient. Detecting both jet and disk gradi¬ 


ents in the same protostellar system remains challeng¬ 
ing, with the few successes to-date leading to no clear 
conclusion. Some studies have reported a match in gra¬ 
dient direction, while others claim no match and others 
still claim no detection of gradients in the jet. Stud¬ 
ies close to the jet footpoint usually require the jet to 
be optically visible, and so Glass H T Tauri stars be¬ 
come the targets. In these cases, matchin g jet gradi¬ 
ents within 100 AU are rep orted for DG Tau (ITesti et al.l 
120021 : IBacciotti et al.l 1200^ and GW Tau (Dougados et 
al. in prep.), but no clear detection of jet rotation is 
reported for HH 30 (iGoffev et al.l[2007l iPetv et al.ll2006l 

iGuilloteau et al.l 1200811 . an^ mis-matched gradients are 

repor ted for RW Aur (|Cabrit et al.l 1200^ iGoffev et al.l 

[20III . 

Persevering with this investigation, our latest study 
examines the jet from RY Tau, a 2 Mq T Tauri star 
of spectral type F8 located at 140 pc. An blueshifted 
atomic jet at position an gle 294° has been detected 
both on large scales i n Ha (ISt-Onge et al.l 1200811 and on 
small scales in [O I (lAgra- Amboage et alJ 1200^ . The 
disk has been reso lved in millimetric continuum emission 
(|Isella et al.ll201(Tf l. These observations constrain the in¬ 
clination of the disk rotation axis to the line of sight 
between 65-75 consistent with the UX Ori behavior. 
We present a high resolution spectro-imaging study with 
GEMINI/NIFS-(-Altair of the RY Tau jet in the [Fe il] 
1.64 fj,ia emission line. We also present a study of the 
disk kinematics conducted with IRAM/PbBI in millimet¬ 
ric CO lines. Together these data sets allow us to search 
for velocity gradients both in the jet and disk, and si¬ 
multaneously perform the two consistency checks for jet 
rotation signatures discussed before. Furthermore, the 
imaging capabilities of NIFS/GEMINI provides an in¬ 
sight into the impact of variability within the jet on the 
derivation of rotation signatures. Our observations are 
described in section [51 while the method of analysis and 
results are detailed in sections |31ll] and |SJ In section [51 we 
discuss the measurements in the context of the magneto¬ 
centrifugal acceleration process, and we summarise our 
conclusions in section |3 

2. OBSERVATIONS 

To observe the outflow from RY Tau (R=9.67), GEM- 
INI/NIFS+Altair was used to obtain AO-corrected (with 
natural guide star) IFU data of the approaching jet 
(PA=294°), with the H-band filter/grating at slit/slice 
position angles both parallel (24°) and anti-parallel 
(204°) to the jet axis (Program ID: GN-2009B-Q-43). 
This configuration positioned the IFU 0.1"x0.043" spax- 
els with the smallest sampling in the direction perpendic¬ 
ular to the jet axis. Additional arc lamp exposures were 
taken during observations to ensure high accuracy wave¬ 
length calibration. Velocity resolution of the instrument 
was 57 km s“^, but for good signal-to-noise we could 
achieve higher velocity precision in emission line veloc¬ 
ity centroids, via Gaussian fitting (section [4]). Spatial 
resolution with AO correction (measured as the FWHM 
of the stellar PSF core) in the direction across the jet 
is 0'.'I2 (17 AU), and along the jet is 0'.'2 (28 AU). To 
retrieve the faint jet emission close to the star (where 
we wish to probe for jet rotation), accurate subtraction 
of the stellar continuum was required. This necessitated 
inclusion of the star within the 3" x 3" field of view, and 
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hence a reduction in the exposure time (to 10 s) with a 
consequent increase in the number of co-added frames, 
in order to achieve a jet signal-to-noise of typically 35. 
The star was offset from the centre of the field so as to 
observe as much of the jet as possible while still ensuring 
coverage of the entire stellar PSF. With this setup, the 
IFU captured the jet out to 1.7" from the star. One third 
of the time was spent on sky exposures, via the standard 
ABA nodding technique. Unfortunately, some of the ob¬ 
servations with the anti-parallel slit show that the field 
of view (FOV) was re-centered on the star, rather than 
offsetting to the jet, causing a loss of information on the 
outer jet emission. Therefore, these observations were 
not included in the co-adds and hence the jet was not 
observed with the same signal-to-noise at each position 
angle. It subsequently transpired that the anti-parallel 
slit data could not be used as it did not reach sufficiently 
high signal-to-noise. 

To observe the RY Tau disk, J=2-l observations 
of the system were performed with the IRAM Plateau 
de Bure Interferometer on 2007 December 18 using the 
extended C configuration with 6 antennas and baselines 
ranging from 24m to 176m (project QOlO). The 1 mm re¬ 
ceivers were tuned single sideband at 230.538 GHz. One 
correlator band of 40 MHz was centered on the J = 
2-1, implying a channel spacing of 78 kHz or 0.1 km s“^. 
The total telescope time, including phase and amplitude 
calibrators, amounted to 24h with 6 antennas. 

3. DATA REDUCTION 

The NIFS data were processed through the standard 
GEMINI calibration pipeline, but using the additional 
arc lamp files to improve wavelength calibration. It was 
not necessary to perform telluric subtraction during the 
pipeline, because subsequent continuum subtraction took 
into account the telluric features. The output IFU cubes 
for a given position angle were then combined, and con¬ 
tinuum subtraction performed as follows. We selected 
a high signal-to-noise reference photospheric spectrum 
free of jet emission. This reference spectrum was then 
scaled down to the level of continuum in a region near 
the emission line of interest, and was subtracted from 
the jet spectrum. The residual emission line profile is 
not only free of photospheric absorption lines but also of 
telluric features. Note, this procedure only fully corrects 
the telluric absorption formed against the stellar contin¬ 
uum. It does not correct for telluric absorption against 
the extended [Fe II] emission. However, we checked that 
this is a small effect (le ss than 10% of the total [F e II] 
emission intensity) (see lAgra-Amboage et al.l (|2011fl for 
further details). In our data, the brightest jet emission 
line was ]Fe II] at 1.644 /im. Continuum subtraction 
was conducted with due care in this region of the spec¬ 
trum, and was sufficiently accurate to allow retrieval of 
jet spectra as close to the jet base as 0'.'2. Observations 
of the sky were processed through the same pipeline as 
the science observations. In this way, we ensure that the 
sky observations carry the same velocity accuracy as the 
science data. We can then profile fit the sky lines to de¬ 
termine the velocity calibration error, which will be the 
same as the science data. 

The Plateau de Bure data were reduced using the 
GILDA^ software. Standard calibration methods, us¬ 
ing close calibrators, were applied. The phase and am¬ 


plitude calibrators 0528-1-134, J0418-I-380 were observed 
along with 3G454.3 as radio-frequency bandpass calibra¬ 
tor. The low quality data were filtered out. In particular, 
data with phase noise rms above 50° were not used in 
subsequent data analysis. J=2-l images were pro¬ 

duced using the GLEAN algorithm. The clean beam size 
is 174500 X I'.'OlOO at PA=41°. Velocities are expressed 
in the LSR frame with a systemic velocity of 7.8 km s“^. 

4. DATA ANALYSIS 

The individual NIFS spectra showed that the [Fe II] 
line was well fitted with a single Gaussian. Emission line 
profile htting was conducted at every position along and 
across the jet. The resulting velocities were corrected by 
-17 km s“^ to the LSR velocity frame. 

The overall error in centroid velocity was determined 
by measuring the velocities of the OH skylines. To do 
this, four of the brightest sky lines were selected (at 
1.60307,1.61283,1.66924 and 1.69035 ^m) based on their 
close proximity to the [Fe II] 1.644 /im jet emission line. 
The velocities were measured in each spaxel, and their 
RMS was calculated separately for each slicing mirror, 
as this is the direction across the jet. Hence, we could 
identify wavelength calibration drifts which may affect 
our measurement of the transverse radial velocity profile 
of the jet. 

Each RMS value comprises a centroid fitting error and 
a velocity calibration error, only the former of which de¬ 
pends on signal-to-noise. That is, the theoretical stan¬ 
dard deviation, a, in estimating the centroid of a G aus- 
sian distribution is obtained by ([Porter et al.ll200Jl : 

FWHM 

“ 2.35 SNR ^ ^ 

where FWHM in this case is measured in the spectral 
direction, and SNR is the signal-to-noise ratio at the 
Gaussian peak. 

For each skyline, the RMS error was plotted against 
the median signal-to-noise, figure [l] and a curve marking 
the theoretical fitting error (according to equation [T|) was 
added for comparison. Binning (x4 along the slitlet) of 
the brightest skyline (snr=10) allowed extrapolation of 
the error to a signal-to-noise comparable to that of the 
jet emission (i.e. SNR=25). Subtracting in quadrature 
the theoretical fitting error (curve) from the total error 
(data points) allows extraction of the velocity calibration 
error, acai, which in our case averages at 1.6 km s“^ 
(RMS). Using this value, we can now calculate the total 
error in centroid velocity for jet emission of any signal- 
to-noise, by simply adding the velocity calibration error 
(i.e. a cal) in quadrature to the fitting error (i.e. the 
output of equation [T|) . 

For example, for a jet signal-to-noise of 25, and given 
the spectral resolution of NIFS (FWHM=57 km s“^), the 
1 a error on the velocity centroid fitting is 1 km s“^ (in 
line with the lower limit on the sky line RMS, figure [T]). 
Adding this in quadrature to the velocity calibration er¬ 
ror of 1.6 km s“^, yields an overall 3 sigma velocity error 
of 5.7 km s“^. Multiplying by a/S, we obtain a 3 tr error 
on differences in Doppler shift of 8 km s“^. Note that 

1 See http://www.iram.fr/IRAMFR/GILDAS for more informa¬ 
tion on GiLDAb software. 
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Table 1 

GEMINI/NIFS+Altair observations of the RY Tan blue-shifted jet. 


Observation Date 

yiitlets FA 

n 

(Jo-added Exposure l imes 

(s) 

'i'otal 'l ime 
(hrs) 

14 Dec 2009 

24 

3060, 1020, 1530, 2040 

2.125 

15 Dec 2009 

204 

4080, 3060 

1.983 

24 Dec 2009 

24 

2040 

0.567 


this is a lower limit, based on a maximum jet signal-to- 
noise of 25-30. For the jet borders, where signal-to-noise 
can drop to 10-15, this error grows to 10-13 km s“^. Pre¬ 
viously published jet rotation indications are typically on 
the order of 1-15 km s“^. So, in spite of our strategy to 
increase velocity accuracy by including wavelength cali¬ 
bration arc lamp exposures during the observations, we 
still find ourselves operating at the very limit of the in¬ 
strument’s ability to detect differences in Doppler shift 
across the jet. 

Finally, note that the calibration error, CTcaZ, includes 
random spaxel to spaxel fluctuations and systematic 
drifts. We examined such drifts and no systematic trends 
in skyline velocity measurements above iacai >5 km 
s“^ were identified across the slitlets/slicing mirrors (i.e. 
along the jet). Along the slitlets (i.e. across the jet), a 
slight trend was identified. However, the gradient did not 
always have the same slope, and the amplitude was low 
2-3 km s“^), hence it remains well within the overall 
three sigma error bars. Furthermore, it is possible that a 
rotation signature could be mimicked by a contribution 
from an uneven slit illumination, due to the strong gra¬ 
dient in brightness within a given slitlet caused by the 
continuum of the star. (Le. when we compute the trans¬ 
verse radial velocity differences, we take into account the 
fact that the jet axis is wiggling such that uneven slit 
illumination effects will introduce a spurious value, even 
if the PSF is exactly symmetric with respect to the axis 
of the detector which it is not). We estimated the am¬ 
plitude of uneven slit illumination effects on transverse 
radial v elocity measurements using the procedure out¬ 
lined in lAgra-Amboage et al.l (120141) . This was found to 
introduce an effect on the order of 0.5 km s“^ which was 
deemed negligible in view of the overall error bars. 

5. RESULTS 

While the NIFS and PDBI datasets presents a wealth 
of information on the jet and disk morphology and kine¬ 
matics, we choose to restrict our current focus to the 
case for jet rotation. Two separate more comprehensive 
studies of jet and disk are to follow (Coffey et al. in 
preparation; Dougados et al. in preparation), with just 
a brief commentary on the jet given below by way of 
introduction. 

5.1. Jet Morphology and Kinematics 

It is interesting to note the general morphology and 
kinematics of the jet, as revealed by the [Fe II] 1.644/im 
emission line. A contour plot of the jet intensity, fig¬ 
ure [H shows peaks are located at positions around -0"2 
to -0"5 and -0"8 to -F.'l, which indicates a clear presence 
of internal jet shocks (or knots) along the flow direction. 
The jet axis is over-plotted as a dashed central curve. 
It was determined by gaussian fitting the intensity pro¬ 
file across the jet image. At each position in the image, 


Velocity Measurement Accuracy 



5 10 15 20 25 

Signal to noise 


Figure 1. RMS velocities plotted for each slicing mirror (slitlet) 
of each skyline. Each skyline is represented by a different symbol. 
The brightest skyline (snr=10) was then binned (x4) along the 
slitlet and the improved RMS was plotted (see points at snr=25). 
The dashed line represents the theoretical Icr error on the Gaussian 
centroid fit, calculated according to equation [T] 


the radial velocity was obtained by Gaussian fitting the 
spectral emission peak where signal-to-noise > 5. Fig¬ 
ure [2] shows how the contours of the jet intensity relate 
to the velocity colour map. Note how peaks in intensity 
contours coincide with local maxima in the on-axis ve¬ 
locity. In the second knot (at -0"8), it is clear that the 
velocity peak slightly precedes the intensity peak, which 
supports the presumption that we are observing emission 
from shock-excited gas. 

Recall that the spatial resolution (measured via the 
stellar PSF) in the direction across the jet is 0"12 (17 
AU), and along the jet is 0'.'2 (28 AU). Meanwhile, we 
measure a jet FWHM of O'.'18 (25 AU), at -0'.'5 (70 AU) 
projected distance from the star. Therefore, the jet width 
is resolved in this region close to the star. Figure [3| shows 
the transverse radial velocity profile, at four sample po¬ 
sitions along the jet. (Figure Sj shows the radial velocity 
profile across the jet at all positions along the jet.) Along 
the jet, we can distinguish two shapes of transverse radial 
velocity profile. A v-shaped profile appears where the jet 
velocity peaks on-axis, at positions -0'.'3, -0'.'8 to -0'.'9 and 
-1'.'2 to -1'.'3. Again, this seems to coincide with intensity 
peaks of the shocked jet material. Between shocks, we 
see a flat profile, at positions -0'.'4 to -O'.'G, -I'.'O to -I'.'l 
and -1'.'4 to -T.'5. 

To add to the complexity, we also see signs of wiggling 
of the jet as it propagates, as indicated by the dashed 
curve in figure [D and the dashed central line of figures [3] 
and m Indeed, the jet wiggle, which becomes apparent 
after -Qt.'b, implies that any kinematic record of the jet 
launch mechanism is probably contaminated beyond this 
point. Nevertheless, in this case, its amplitude is small 
and we take it into account when considering transverse 
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Figure 2. Radial velocity map of the RY Tau jet in [Fe II] 1.644 
/im emission. The map is overlaid with intensity contours, and the 
position of the jet axis (dashed line). The approaching jet PA = 
294 °. 

kinematic signatures which may represent jet rotation. 
Although the wiggle seems to cause the jet to veer to the 
side at -1"25, the jet is in fact opening out into a v-shape, 
the left arm of which falls below the 5 sigma level. 

Our data show that radial velocity profiles can be dom¬ 
inated by shock kinematics and jet wiggle, rather than 
any record of the jet launch mechanism, once the jet has 
propagated beyond about -C/.'b. This highlights the need 
to examine kinematics as close to the source as possible, 
in order to reduce the impact of these effects. 

5.2. Jet Rotation 

If the jet is rotating, we should observe a difference in 
the Doppler shift between one side of the jet axis and the 
other. Examining radial velocity measurements across 
the jet, we do not see any clear sign of jet rotation at the 
3(7 level. In the right-hand panels of figure |31 we compare 
the radial velocities at symmetric positions about the jet 
axis (taking account of the jet wiggle). We do not find 
any clear sign of systematic velocity asymmetries across 
the jet at the 3a level, at any distance from the star. 
In the region close to the star -0"3 to -0"6), we do 
see a consistent slope in radial velocity profile over the 
3-5 spaxals across the jet. Figure [3] shows that, at a 
distance of -0"5 from the star, the peak signal-to-noise 
on the jet axis is 30. Interpolating this plot allows us 
to more precisely examine matching positions either side 
of the jet axis. Table [2] give the results at -0"5 from the 
source where the velocity is largest. We find differences in 
Doppler shift of 7±3 km s“^ (Icr) at a radius of 0'.'07 (i.e. 
10 AU) from the jet axis where the signal-to-noise drops 
to 19. Further from the axis, at O'.'l, the signal-to-noise 


Table 2 

Differences in Doppler Shift across the RY Tau jet at 0(^5 from 
star, based on interpolation of data in figure [3] Errors are at the 
1(7 level. 


Radius trom axis 

SNR 

Vi 

V2 

V2 - Vi 

(arcsec) 


(km s“^) 

(km s“^) 

(km s“^) 

007 

19 

-71±2 

-78±2 

-7±3 

0.1 

11 

-72±4 

-77±4 

-5±6 


drops even more (to 11) and the difference is 4±6 km s“^ 
(Icr). Hence, although it is possible that a signature of 
jet rotation in present in the innermost jet region, it is at 
most at the 2cr level and thus not a conclusive detection. 
Unfortunately, the jet signal-to-noise for the anti-parallel 
slit was not high enough to verify this trend. The slope 
persists until -0'.'6, just before the first jet knot, but not 
beyond. 

Jet radial velocities can be used to find the jet incli¬ 
nation angle, and hence the jet toroidal velocity. Radial 
velocities are in the range -60 to -110 km s“^. We note 
that the radial velocity is highest (^ -90 to -108 km s“^) 
in the region closest to the source (^ -073), and at var¬ 
ious points along the jet axis. Taking a typical value 
of -75 km s~^, and a proper motion of ^ 165 km s“^ 
(|St-Onge et al.l[2008f) . we calculate an inclination angle 
with respect to the line of sight of 65.5°. Note that this 
rough estimate of the jet inclination from our observed 
radial velocities and the proper motions reported on large 
scales is compati ble with the disk inclination on 100 AU 
scales of 70±5° (jlsella et al.ll2010fl . We use this to find 
the jet toroidal velocity (section!^. 

5.3. Disk Rotation 

Plateau de Bure observations of the RY Tau disk re¬ 
veal detections of ^^CO(2-l) emission from vlsr = l-fO 
to 11.55 km s“^, i.e. covering a velocity range of 10.3 
km s“^. Contour plots of the CO red and blue wing 
emissions, figure [5l reveal a clear velocity gradient along 
PA=24°, i.e. perpendicular to the jet axis and con sistent 
with the disk PA derived by llsella et al.l (|2010D . The 
sense of the velocity gradient (blue-shifted emission to 
the north and red-shifted to the south) implies that the 
disk spin vector points along the axis of the approaching 
jet. A contour plot of the CO position-velocity diagram 
along PA=24° shows that the velocity held is consistent 
with Keplerian rotation around a 2 Mq star, hgure O 
The value of the systemic velocity suggested by the CO 
kinematics is estimated as Vgys = 6.75 ± 0.2 km s“^. We 
use an inclination o f the disk rotatio n axis to the line 
of sight of 70 ± 5 ° (jlsella et al.ll20inD . A more detailed 
study of the disk, via more complex modelling, is beyond 
the scope of this paper and is to be presented separately 
(Dougados et al. in preparation). 

6. DISCUSSION 

In this section, we examine the implications for jet 
launching, under a steady state assumption, of our de¬ 
rived upper limits on radial velocity differences across 
the jet. 

6.1. Upper limits on jet toroidal velocity and angular 
momentum at 70 AU above the disk-plane 
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Figure 3. Left panels: Transverse velocity profiles along the jet at various positions: -0.5”, -0.8”, -1.1” and -1.3”. The position of the jet 
axis is indicated by the dashed line. Right panels: comparison of centroid velocities on either side of the jet axis. Error bars are at the Itr 
level. 


While we do not detect any consistent rotation signa¬ 
ture above the 3 a level in any of our transverse cuts 
of the jet, we can use our 3a upper limit on (v 2 — vi) 
derived in section [Sj to determine upper limits on any 
possible toroidal velocity that may be present. 

In doing so, we note that an important consideration 
when attempting to constrain models for jet launching 
with signatures of jet toroidal velocities is that measured 
velocity differences across the jet always underestimate 
the true rotation profile, especially at small transverse 
distances from the jet axis. This is due to integration 
of the flow along the line of s ight, combined with beam 
smearing (jPesenti et al.l [200^ . Correction for this de¬ 
pends on the spatial resolution of the observations, the 
distance from the jet axis where differences are measured, 
the jet velocity interval considered, and the jet ionisa- 
ti on profile, among ot her parameters. From the study 
of iPesenti et al.l (j2004li . we estimate that beam smearing 
effects can underestimate the toroidal velocity by up to a 
factor 2 at r=0"07, and a factor 1.2 at r=(y.'l. In order to 
minimise the beam smearing effects, we thus take the 3 
a value at r=0'.'l to estimate the upper limit on {v 2 —vi) 
to find an upper limit on the toroidal velocity. 

For a jet inclination angle of i =70° to the line of sight, 
the 3 a value at r=0'.'l of 18 km s“^ translates to an upper 


limit on the toroidal velocity of = {v 2 — vi) / {2sin(i)) 
< 10 km s“^ at a radius from the jet axis of r = 14 AU 
and a distance above the disk-plane of z = 70 AU. This 
implies a maximum value for specific angular momentum 
of ru0=14O AU km s“^. 

6.2. Implications for jet launching under a steady state, 
axisymmetric assumption 

In examining the implications of our toroidal velocity 
upper limit on jet launching models, we first consider 
that the steady state assumption is valid for the RY Tau 
jet at 70 AU from the source. Although shocks and wig¬ 
gles are clearly present along the jet, they may only repre¬ 
sent small perturbations in an underlying steady launch¬ 
ing process. I.e. since the ambient medium has already 
been cleared by the jet in this evolved T Tauri system, 
we can be sure that the shocks seen in [Fe II] are internal 
shocks where high-velocity jet material meets previously 
ejected lower velocity jet material. Hence, the observed 
radial velocities represent the jet radial velocity in the 
post-shock regions. Indeed, observed proper motions and 
analyses of line ratios in protostellar jets indicate shock 
speeds less t han 25 % of the flow velocity on scales of a 
few 100 AU (|Lavallev-Fouauet et al.l[2000li . Such shock 
amplitudes are consistent with the observed variation of 
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Figure 5. ^^CO(2-l) channel maps from Plateau de Bure obser¬ 
vations, illustrating the disk rotation sense. Blue and red contours 
show emission integrated from LSR velocities of -2.8 to -4.4 

km s“^, and 9.5 to 11.1 km s“^, respectively. Contours start at 3 
<7. The beam size is shown as an insert. A clear velocity gradient 
is observed at PA=24°, perpendicular to the black arrow which 
represents the approaching jet (PA= 294°). 



Figure 6. ^^CO(2-l) position-velocity diagram along the disk PA, 
from Plateau de Bure observations. The two green curves show 
predictions for Keplerian rotation around a 1 and 2 Mq star. The 
overlap of the contour peaks with the curves illustrates that the 
disk is in Keplerian otation . A disk inclination of 70° was used, 
based on llsella et al.l II2Q10I) . The dashed line marks an estimate 
of the systemic velocity of 6.75 km s“^. 

radial velocities in the inner RY Tan jet, i.e. {v 2 — vi) 
= 20 km compared to the mean of 85 km s“^. This 
also suggests that the observed flow velocities should be 
close to the mean ejection velocities from the disk. 

Toroidal velocity is a valuable parameter in differenti¬ 
ating between com peting steady-state MHD launching 
mechanisms (e.g. iByciotti et al.l l2002t iFerreira et al'l 
120061 : I Anderson et al.ll2003D . Under the assumption of 
steady MHD ejection, a relation can be found between 
the specific angular momentum(ru 0 ) and the poloidal ve¬ 
locity (vp) for a given jet streamline. This relationship 
depends on one parameter: in the case of disk winds, the 
magnetic lever arm, A, or launching radius , tq, such that 
(|Ferreira et al.ll200a lAnderson et a.l.ll2003ll : 

^ (2) 



Vp / (MVMq)'^ kms-‘ 

Figure 7. Comparison of launching mechanisms as a function of 
the relationship between jet p o loidal and specific angular momen¬ 
tum (similar to IFerreira et aLl ||2006I) V RY Tau is plotted accord¬ 
ing to our measurements. The RY Tau coordinates imply an upper 
limit on the jet launch radius in the disk of ro << 0.45 AU, under 
the assumption that the jet is still in steady-state at 0(^5 (70 AU) 
from the star and a value of < 5. Recall that the range in Vp 
comes from the range for the jet inclination angle of 70 ±5°. 

2rv^no = vl + (3) 

where Hq is the angular velocity at a radius in the disk- 
plane of ro, and is the effective lever arm at a given 
position from the driving source. Note that the value of 
will remain lower than the true magnetic lever arm 
value, A, until all angular momentum has been trans¬ 
ferred to the matter. As noted in citealpFerreira2006, 
this complete transfer may be reached only quite far from 
the star. For example, they show a solution with A = 13 
which has A^ = 8 at as distance from the source of 50 
AU. From these two equations, and assuming no pres¬ 
sure effects (/3=0), we see that simultaneous knowledge 
of both Vp and allows us to constrain both vq and 
Xfj). We stress that these derived relations are very gen¬ 
eral, and are valid for all types of super-Alfvenic, station¬ 
ary, axisymmetric, self-collimated MHD jets, since they 
only express the conservation of total angular momen- 
tum and ener g y alon g a given streamlines. We refer to 
IFerreira et al.l (1200611 for details of how to derive these 
equations. 

Figure[7]illustrates the parameter space of the ver¬ 
sus Vp plane. Plot contours highlight values of A^ and 
ro. Only predictions for extended disk wind models are 
shown. The X-wind is a li miting case of these solution s 
with ro=0.05 AU and A=3 (lAgra-Amboage et al.ll2009tl . 
Stellar wind predict ions fall below the lowest vq curve 
(jFerreira et al.ll2(i0^ . 

A smaller toroidal velocity at a given radial distance 
from the jet axis implies smaller specific angular momen¬ 
tum, rv^. Figure [7] shows that for a given Vp, this then 
implies a smaller launching radius, vq. Hence, we can use 
our derived upper limit on the specific angular momen- 
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turn to derive an upper limit on the launching radius of 
the [Fe II]-emitting jet from RY Tau. 

Using our velocity measuremen ts, we plot RY Tau on 
the graph bv iFerreira et al.l (|2n06[l . To do this, we calcu¬ 
late the magnitude of the jet poloidal velocity, (from the 
radial velocities of -71 to -78 km s“^ from table[2]and jet 
inclination angle of i=70±5°) to be -168 to -301 km s“^, 
and take the stellar mass of M^^/Mq = 2. The result¬ 
ing coordinate positioning indicates an upper limit on 
the launch radius in the disk, ro< 0.45 AU, and a small 
magnetic lever arm parameter, A=(r^/ro)^< 7. Note 
that our constraint on rp is only applicable for the part 
of the flow with radial velocities larger than the value 
at the radius from the jet axis of r=0'.'l and at 70 AU 
from the source. There are clearly data points with lower 
radial velocities in the outer borders of the jet, and this 
jet material could be originating from larger radii than 
rg <0.45 AU. However, in our observations, these data 
are too faint to analyse. 

We briefly consider the implications of mass fluxes for 
the launching radius in the RY Tau system. In a steady 
MHD disk wind that carries away all of the accretion 
angular momentum flux from a Keplerian disk, there is 
a single inverse relation which has to hold between the 
magnetic lever arm parameter, A, and the ratio of ejected 
to accreted mass-flux. We can then check that the de¬ 
rived \,p and ro paramete rs are compatible wi th the ob¬ 
served jet mass loss rate (|Ferreira et al.l [200^ equation 
17): 


A/j 1 ^ out 

Ma 4(A 1) Tin 


( 4 ) 


where Mj and Ma are the mass ejection and accretion 
rates respectively, and Tin and Tout are the inner and 
outer ejection radii on the disk plane. It is not a degen¬ 
eracy of oversimplified self-similar models, but a basic 
requirement coming from en ergy and angular m omen- 
tum conservation (see e.g. iPudritz et al.l 1200711 . An 
estimate of the mass flux in the RY Tau blue-shifted 
atomic jet is between 0.16 and 2.6 x 10 “® Mq yr“^ 
l|Agra-Amboage et al.l 1201 Hi . while accretion rate esti¬ 
mates range between 6.4-9 x 10“® Mq yr“^. This gives 
the ratio Mj/Macc = 0.02-0.4. If we take our constraint 
on rout = ro < 0.45 AU, and assuming =0.1 AU 
(i.e. close to co-rotation), then the (one-sided) mass flow 
ratio in RY Tau implies A < 1.9-19.8. This is consistent 
with our derived constraint A < 7. 


6.3. Implications for jet launching without the steady 
state, axisymmetric assumption 

We stress that our upper limits on rg and A^ are only 
valid if the jet is in steady state. Otherwise, the MHD 
invariants used to derive the formulae above would no 
longer be applicable. Indeed, the steady-state assump¬ 
tion may be uncertain, even at 075 from the source, since 
our intensity map (fig [2]) shows an emission knot as close 
as 0735 from the star which is associated with a velocity 
peak at -110 km s“^, indicating the occurrence of shocks 
very close to the source. Similarly, the detection of small 
jet wiggles beyond 074 also makes the assumption of ax- 
isymmetry uncertain at this position. 

Recent numerical works have investigated the effect of 
shocks and jet propagation on rotation signatures in disk 


winds iFendtl (120111) show that magneto-hydrodynamic 
shocks occurring in an initially non-rotating jet with he¬ 
lical magnetic field configuration can transfer magnetic 
angular momentum (magnetic stress) into k inetic angu¬ 
lar m omentum (jet rotation). Conversely, iSaiitv et al.l 
(I2012D demonstrate that layers of a rotating disk wind 
may be temporarily counter-rotating due to initial varia¬ 
tions in the flow velocity, without contradicting the MHD 
jet launching mechanism. Similarly, Staff et al. (2014, 
2015) report 3D MHD simulations of disk wind propa¬ 
gation which show that kink instabilities can result in a 
signature of rotation in the outer jet layers in a direc¬ 
tion opposite to that of the disk. If indeed time depen¬ 
dent effects are critical in structuring jets even as close 
to the star as investigated here (i.e. 70 AU), then our 
derived upper limit on angular momentum does not pro¬ 
vide a reliable constraint on the launching radius and 
magnetic lever arm for RY Tau. Smaller spatial scales 
(^ 10 AU) should be investigated, where either the jet 
is steady and axisymmetric or where direct comparison 
with time-dependant 3D MHD numerical simulations in¬ 
cluding disk feedback becomes possible. 

7. CONCLUSIONS 

High resolution observations of jets close to their base 
are critical in constraining the jet ejection mechanism. 
Here, we present results for RY Tau, via an IFU dataset 
of high resolution both spatially and spectrally. We also 
present Plateau de Bure observations of the disk, which 
reveal a velocity field consistent with disk rotation, and 
we attempt to determine if we can observe a transfer 
of angular momentum from disk to jet. A radial velocity 
analysis of the jet revealed that there is no systematic dif¬ 
ference across the jet at the 3 a level, i.e. no indication of 
jet rotation. Nevertheless, we can use the 3 a error bars 
to obtain a useful upper limit on the current jet toroidal 
velocity. For RY Tau, at 075 (70 AU) from star and 07l 
(14 AU) from axis, we find an upper limit of 10 km s“^ 
on the toroidal velocity. If we assume the jet is still in a 
steady-state at this distance from the star, the toroidal 
velocity would in turn provide an upper limit on the jet 
ejection radius in the disk plane of rg < 0.45 AU, and a 
magnetic lever arm parameter of < 5. (Recall that our 
constraint on rg does not hold for jet material of lower 
radial velocity in this region which may originate from a 
larger launching radius.) These results cannot differen¬ 
tiate between competing models but, strictly under the 
assumption of steady-state, the measurements can place 
important constraints on steady-state MHD models of jet 
launching. However, the steady-state assumption is chal¬ 
lenged by our findings of multiple internal shocks start¬ 
ing as close as 073 (50 AU) from the source, and a jet 
wiggle which develops further out. In this case, no firm 
constraint can be placed on the jet launch point since 
most of its angular momentum could be stored in mag¬ 
netic form, rather than in rotation of matter. Although 
we do not reach the velocity accuracy anticipated, the 
variations within the IFU dataset make it clear that a 
three-dimensional datacube is critical to interpretations 
of velocity asymmetries across T Tauri jets. 
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